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Hierarchy of gaps and magnetic minibands in graphene 
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We determine bands and gaps in graphene subjected to the magnetic field of Abrikosov lattice of 
vortices in the underlying superconducting film. The spectrum features one non-dispersive magnetic 
miniband at zero energy, separated by the largest gaps in the miniband spectrum from a pair of 
minibands resembling slightly broadened first Landau level in graphene, suggesting the persistence 
of = ±2 and ±6 quantum Hall effect states. Also, we identify occasional merging point of magnetic 
minibands with a Dirac-type dispersion at the miniband edges. 

PACS numbers: 73.22.Pr, 73.21.Cd,73.43.- 


Studies of superlattices in two-dimensional (2D) elec¬ 
tron systems have, recently, been boosted by the devel¬ 
opment of van der Waals heterostructures of graphene 
with hexagonal boron nitride (hBN). In such systems 
the superlattice effects, observed in STM spectra m, 
magneto-transport chracteristics [IHg and quantum ca¬ 
pacitance [7], are produced by a periodic moire pattern, 
with the period a determined by slight incommensurabil¬ 
ity and misalignment between graphene and hBN crys¬ 
tals [DEIS! and reflect the formation of superlattice 
minibands for graphene’s Dirac electrons EEIIS]. To a 
large extent, the possibility to observe the superlattice ef¬ 
fects in graphene-hBN heterostructures owes to the high 
mobility of electrons in such systems, where graphene is 
encapsulated between hBN sheets both protecting from 
contamination and permitting to vary electrons’ density 
over a broad range using electrostatic gates. When sub¬ 
jected to a strong external magnetic field, the superlat¬ 
tice leads to the formation of a ’Hofstadter butterfly’, a 
sparse spectrum of minibands [inUT^ formed at magnetic 
field values corresponding to the magnetic flux, ^ — ^(j)Q 

(through the area S = 'JZcP' {2 of the superlattice unit 
cell) commensurate with the flux quantum, (jjQ = hje. 

Here, we consider a magnetic superlattice [TMIS] that 
can be realised in a ballistic hBN-graphene-hBN stack by 
placing it over a high-7Lc2 superconductor film {e.g., Nb, 
W, or MoRe alloy). In such a system, where no alignment 
control of graphene and hBN lattices is required, long- 
range periodic structure is caused by the Abrikosov lat¬ 
tice of vortices muD] formed in a superconductor sub¬ 
jected to an external magnetic field H < iJc 2 , sketched 
in the inset in Fig. In contrast to the earlier theories 
developed for spatially alternating magnetic fields with a 
zero average [HHH] , the Abrikosov lattice produces mag¬ 
netic induction with spatial average, B = (fto/{'\/ia?), 
linked to the magnetic lattice period, a. As each vor¬ 
tex carries the flux /i/2e, the vortex lattice realises the 
simplest fundamental fraction, | = i, in the Brown-Zak 
commensurability condition for magnetic field flux in a 
2D periodic system mm- 

Figure [l] shows the hierarchy of bands and gaps in 
the corresponding spectrum of Dirac electrons calculated 



FIG. 1: Spectrum of Dirac electrons in graphene in a magnetic 
field of Abrikosov vortex lattice, with one degenerate band 
at F = 0. Energy is scaled as vy/2heB = ja, and 

magnetic field as A-k}?B j<f>Q. Inset: 2nd and 3rd minibands 
dispersion over the folded magnetic Brillouin minizone near 
their touching condition. 


in this work and plotted as a function of 47rA^H/^o- 
The latter parameter characterises the the ratio between 

the lattice period a = -sj4>o/\/iB and the penetration 
depth A in a superconductor. The spectrum in Fig{^ 
features one degenerate magnetic miniband which pre¬ 
cisely coinsides with the zero-energy m = 0 Landau level 
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(LL) peculiar for Dirac electrons, and two other low- 
energy bands which resemble slightly broadened m = ±1 
LLs. At the same time, the higher-energy minibands 
(traceable at inX^B/ipo > 1 to LLs with \m\ > 1 at 
^ ~ |^^^\/2|m.|fiei3, m S Z) are strongly broadened, 
and they overlap on the energy scale when a ^ A, form¬ 
ing Dirac-type features at conjoint miniband edges. 

Distribution of magnetic field of Abrikosov’s lattice 
(for isotropic superconductors, the vortex lattice is 
hexagonal) is given by [THl HO] : 

iL(r) = ^iL„(r-R,), 

where Ri picks location of each individual vortex with 
the field profile given by the modified Bessel function 
of imaginary argument, with Ko{x <?; 1 ) « In A and 
Ko{x S> 1) « convenience, we use Fourier 

representation for the periodic field of vortex lattice [ 20 ] 
and the corresponding vector potential. 


To find the magnetic miniband spectrum correspond¬ 
ing to Hamiltonian we use [12] the basis of Bloch 
states (with s = —Nj2, ■ ■ ■ ,JV/2, and t = 0, 1 ), 
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and a non-orthogonal coordinate system adjusted to the 
hexagonal symmetry of the vortex superlattice. 


Here, Mm are Hermite polynomials, and the sign of m 
identifies the conduction {m > 0 ) and valence (to < 0 ) 
band states. 

The states in the basis set |5f‘(fc)) transform according 
to the irreducible representations of the symmetry group 
A4g of the vortex lattice field, which includes Cg rotations 
and magnetic translations ^ = {0x = -X' = 

n'laxi +n 2 ax 2 } C A4g. In contrast to usual translations, 
magnetic translations do not commute with each other, 
0 ai 0 a 2 = —©asQai, howevei', the group contains 
an abelian subgroup. 


Snina = 5^ X (’^•1*1 + ?T-2®2), 


w„i„2 = ynf+ n|+nin2, g = 47r/(V3a) 

1 . V3 . . 1 . V3 . 

xi = + X 2 = --x+—y, 


The spectrum of electrons in K and K’ valleys of 
graphene is determined by the Hamiltonian 
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where A = hxi{—xi -\- 2x2)/{ia^e), xi = x + 

X 2 = —X + and v « 10®cm/s is Dirac velocity in 
graphene. This Hamiltonian acts in the space of two- 
component wave functions describing electrons’ ampli¬ 
tudes on A and B sublattices of the honeycomb lattice of 
carbons, with the basis choice ['I'(A),'I'(_B)] in valley K 
and ['I'(i3), —'I'(A)] in valley K’: this choice provides us 
with the same form of the Hamiltonian in both valleys. 


g' = {Qji = R = 2 n[axi + 271^X2}, 

which is formed by translations on a superlattice with 
isotropically doubled period and a unit cell area 45. Be¬ 
cause of this, it is possible to classify the states on ma- 
gentic superlattice using the wave vector q taken over 
the folded Brillouin zone (BZ) with the area ^ 5 ^, four 
times smaller than the BZ area of the gometrical vortex 
lattice. Each of these folded states is two-fold degenerate 
[101 I12j . which is prescribed by the anti-commutation, 
0 ai 0 a 2 — 0 ot 2 ®^ii 5 of the operators of elementary 

translations. By analysing characters of the magnetic 
translation group A4§, one can find that the latter fea¬ 
tures six different 2 -dimensional irreducible representa¬ 
tions related to the states with the wave vector <7 = 0 in 
the centre of magnetic BZ. In practice, these six types of 
irredicible representations can be constructed using lin¬ 
ear combinations of Bloch functions \ with 

different M >0 but fixed A^ = 0,1, 2,3,4, 5. 

Using the basis of Bloch functions in Eq. (§ , Hamil¬ 
tonian Eq. ([^, can be represented in the form of the 
Heisenberg matrix. 
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where z„ = + aa). 

In Figure we show bands (shaded) and gaps (white 
intervals) in the spectrum obtained by numerical diago- 
nalisation of Heisenberg matrix In this calculation 
we used 80 LLs to guaratee the convergence of the en¬ 
ergies in the lowest 20 bands on the conduction and va¬ 
lence band sides, and we included all points gmn^ in the 
reciprocal space such that Ififumal ^ 647r/(73a). The 
calculated energies are scaled with vV2HeB = ^^Hv/a 
(energy of |m| = 1 LL an homogeneous field B), and they 
are shown for various values of ATrX^BjcjiQ, parameter 
we use to characterise magnetic field distribution across 
the unit cell. The band diagram in Fig. [l] is electron- 
hole symmetric, and the hierarchy of bands and gaps in 
this spectrum is universal, as it can be applied to Dirac 
electrons in the field of a vortex lattice in a film of any 
isotropic type-II superconductor. To mention, due to the 
large demagnetisation factor of a film, vortices enter the 
film at the field much lower than the nominal first critical 
field of a bulk superconductor, hence, justify the regime 
of AirX^BfcjiQ 1. We have set the upper side of the 
interval of inX^B/cjjo < 2 shown in Fig. [^at the lowest 
possible limit for the second critical field Hc 2 , knowing 
that at higher values of Att\'^B/4)q the minibands would 


converge further towards the LL spectrum in a homoge¬ 
nous magnetic field. 

The spectrum in Fig. features a dispersionless zero- 
energy band that exactly coincides with m = 0 LL in 
graphene, both by its energy, large gaps separating it 
from the rest of the spectrum, and capacity to accom¬ 
modate electrons. Also, the first magnetic minibands, 
both on conduction and valence band side, are only 
slightly narrowed and follow almost exactly the energy 
of TO = ±1 LL. All the other magnetic minibands pro¬ 
gressively broaden upon the increase of parameter a, 
and some touch at certain values of a. Such band de¬ 
generacies occur due to occasional non-avoided crossings 
of q = 0 energy levels at the mini-BZ centre (y-point). 
These occasional crossings are allowed by symmetry, be¬ 
cause the energy levels at the mini-BZ centre belong 
to one of six different irreducible representations, so that 
the closest levels do not necessarily ’repel’ each other on 
the energy axis. 

As different irredicible representations are built by 
mixing parent LL states which are, at least, \5m\ = 6 
apart, such crossings can be identified by analysing only 
the diagonal entries in the Heisenberg matrix ([^, 
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where iFi is Hummer’s confluent hypergeometric func¬ 
tion. Hence, we estimate that, within the interval a < 2 
such crossings whould appear at 02,3 = 0.98 for bands 2 
and 3, at 03,7 = 0.90 for bands 6 and 7, and 03445 = 1.03 
for bands 14 and 15, and these values are close to the 
band touching points found in the exact diagonalisation 
of Eq. . At a finite q, such separation of the spectrum 
into six independent groups is no more possible. 

The separation vs mixing of the groups of levels can be 
followed for the states at q tt/c near 7 -point, where 
off-diagonal matrix elements between the closest energy 
states can evaluated analytically [55], leading to a 2 x 2 
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mials (m 7^ 0,±1). This matrix catches the Dirac-type 
edges of touching bands at a = am,m+i, confirmed by 
the numerically calculated dispersions plotted over the 
entire magnetic BZ in the inset in Fig. Note that 
velocity of all these ’second generation’ Dirac electrons 
appears to be almost the same, Vm,m+i ~ 0.75?;, and, 
due to the spin and valley degeneracy which has to be 
factored additionally to the two-fold degeneracy of states 
in magnetic minibands at (/> = ^4>o, each of the calcu¬ 
lated orbital states in the folded magnetic BZ is 8-fold 
degenerate. 



S S 


FIG. 2: Distribution function P{s) of the normalised level 
spacings in the miniband energies at 7 and p points of the 
Broillouin minizone. Energy levels were taken from more than 
30 states sampled (E = 0 and the 1st band excluded) for 10 
values of 0.5 < a < 1.5, with a step of &a = 0.1. 

The mixing vs separation of subsets of states in the 
magnetic minibands can also be traced using the distri¬ 
bution function P{s) of the normalised level spacings. 


Fig. [^for two given points in the Brillouin minizone. For 
the /r-point, it shows strong ’level repulsion’ character¬ 
istic for the unitary symmetry class of random matrix 
theory [5D]. For the 7-point, levels can appear close to 
each other, as happens in other periodic systems 
I34j where high lattice symmetry splits the spectrum into 
subsets of states corresponding to different irreducible 
representations of the lattice symmetry group which can 
appear arbitrarily close to each other [HHHM] . 

The features we found for the Dirac electrons in 
hBN-encapsulated graphene placed over a hexagonal 
Abrikosov vortex lattice suggest that the quantum Hall 
effect at the filling factors ly = ±2 and, to some extent, 
v = ±6 would remain a robust feature in the transport 
and capacitance measurements, whereas Shoubnikov - 
de Haas oscillations at higher filling factors would be 
strongly suppressed. This property of Dirac electrons in 
graphene seems to be generic for a broad range of mag¬ 
netic field distributions. To stress this point, in Fig. 
we show the calculated spetrum of Dirac electrons mov¬ 
ing in the a square vortex lattice, which has all the same 
characteristic features as the spectrum corresponding to 
the hexagonal vortex lattice. 
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